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ABSTRACT 

The formation of large scale structure is independent of the nature of the cold dark 
matter (CDM), however the fate of very small scale inhomogeneities depends on the 
micro-physics of the CDM particles. We investigate the matter power spectrum for 
scales that enter the Hubble radius well before matter-radiation equality, and follow 
its evolution until the time when the first inhomogeneities become non-linear. Our 
focus lies on weakly interacting massive particles (WIMPs), and as a concrete exam- 
ple we analyze the case when the lightest supersymmetric particle is a bino. We show 
that collisional damping and free-streaming of WIMPs lead to a matter power spec- 
trum with a sharp cut-off at about lO^^Af© and a maximum close to that cut-off. We 
also calculate the transfer function for the growth of the inhomogeneities in the lin- 
ear regime. These three effects (collisional damping, free-streaming and gravitational 
growth) are combined to provide a WMAP normalized primordial CDM power spec- 
trum, which could serve as an input for high resolution CDM simulations. The smallest 
inhomogeneities typically enter the non-linear regime at a redshift of about 60. 
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1 INTRODUCTION 

Analysis of the anisotropies in the cosmic microwave back- 
ground (CMB) radiation (Spergel et al. 2003) finds that the 
relative matter density fim = 0.29 ± 0.07 is significantly 
larger than the relative baryon density Q,h = 0.047 ± 0.006. 
These results are consistent with the observed abundances 
of light elements and primordial nucleosynthesis (see e.g. 
Tytler et al. (2000)) and the power spectrum found from 
galaxy red-shift surveys (Percival et al. 2001), and indi- 
cate that the Universe contains a significant amount of non- 
baryonic cold dark matter (CDM). 

Weakly interacting massive particles (WIMPs) are at- 
tractive CDM candidates, since a stable relic from the elec- 
troweak scale generically has an interesting present day 
density, ficdm ~ 0(1) (Dimopoulous 1990). The argument 
goes as follows: annihilation processes cease and the WIMP 
number density n becomes fixed when the annihilation 
rate Fann drops below the expansion rate H (usually re- 
ferred to as chemical decoupling or freeze-out). The tem- 
perature of chemical decoupling is thus defined by Fann = 
{o'!!,nnv)n{Tcd) ~ H{Tcd), where (Jann denotes the annihila- 
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tion cross section and v the relative velocity of the annihi- 
lating particles. For a typical weakly interacting particle one 
finds Ted ~ m/25, m being the WIMP mass. The present 
day relative WIMP density can be estimated as 



mncd(acd/ao)^ 
(SHoVSttG) 



0.2 



(m/rcd)/25 

(crannW)/l pb ' 



(1) 



where a denotes the scale factor and we have set ^ = 0.7 
and assumed 90 relativistic degrees of freedom at chemical 
decoupling. A cross section of about 1 pb (= lO'^^m^) is 
typical for the annihilation of supersymmetric WIMPs since 
1 pb ~ a^/(100 GeV)^ where a 1/100 is the electroweak 
coupling and supersymmetric particles typically have mass 
of order the electroweak scale (100 GeV). 

In supersymmetry every standard model particle has an 
supersymmetric partner and in most models there is a con- 
served quantum number (R-parity), which makes the light- 
est supersymmetric particle stable. Supersymmetry models 
have a large number of free parameters, however in most 
models the lightest supersymmetric particle is the lightest 
neutralino (which is a mix of the supersymmetric partners of 
the photon, the Z and the Higgs bosons; see e.g. Jungman, 
Kamionkowski & Griest 1996), furthermore in large regions 
of parameter space the lightest neutralino is mainly the bino 
(Roszkowski 1991). Accelerator searches place a lower limit 
on the neutralino mass of m > 37 GeV (see e.g. Hagiwara 
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et al. 2002), while the WMAP measurement of the matter 
density leads to an upper limit of m < 500 GeV (Ellis et al. 
2003). 

In CDM cosmologies structure forms hierarchically; 
galaxy haloes form from the merger and accretion of sub- 
haloes (which themselves formed from smaller subhaloes). 
The internal structure of galaxy haloes is determined by the 
dynamical processes which act on the accreted components. 
Dynamical friction causes subhaloes with mass M > 10^ Mq 
to spiral toward the centre of a Milky Way mass parent halo 
within a Hubble time, while the tidal field of the parent 
halo can strip material from a subhalo. Numerically sim- 
ulated galaxy haloes contain large numbers of subhaloes 
that have not been destroyed (Klypin et al. 1999; Moore 
et al. 1999) and the subhalo mass function varies roughly 
as N{M) oc M~^ ® down to the resolution limit of the sim- 
ulations, M ~ W^Mq (Stoehr et. al. 2003). Furthermore 
the anomalous flux ratios of multiply imaged QSOs (Mao & 
Schneider 1998; Schechter 2003) may be due to millilensing 
by subhaloes, providing tentative observational support for 
the existence of substructure in galaxy haloes. 

Substructure on scales far smaller than those resolved 
by numerical simulations has potentially significant conse- 
quences for WIMP direct (Drukier, Freese & Spergel 1986) 
and indirect (see e.g. L. Bergstrom 2000) detection. WIMP 
direct detection (using terrestrial detectors) probes the dark 
matter distribution on sub-milli-pc scales (Silk & Stebbins 
1993; Moore et al. 2001; Green 2002, 2003). WIMP indi- 
rect detection via WIMP annihilation products (gamma- 
rays, antiprotons and neutrinos) is most sensitive to the 
highest density regions of the Milky Way (e.g. Silk & Steb- 
bins 1993; Bergstrom et al. 1999, 2001; Calcarioo-Roldan 
& Moore 2000) and clumping will also enhance the extra- 
galactic gamma-ray signal (Ullio et al. 2002; Taylor & Silk 
2003). 

In this letter we present the matter power spectrum 
close to the end of the linear regime, for scales that enter 
the horizon well before matter-radiation equality. Three ef- 
fects are important: after chemical decoupling the primordial 
CDM density perturbations are collisionally damped, due to 
elastic interactions with other species, and then after kinetic 
decoupling free-streaming leads to further damping. CDM 
inhomogeneities on physical scales larger than the damping 
scale grow logarithmically during the radiation dominated 
epoch and then roughly proportional to the scale factor dur- 
ing the matter dominated epoch. 

The small scale CDM power spectra for the density, 
velocity and potential perturbations are an essential in- 
put for attempts to analytically estimate the size and 
mass of the very first gravitationally bound objects. Ul- 
timately, very high resolution CDM simulations will be 
necessary to make detailed predictions for the fate of the 
very first objects, however the very difi^erent character- 
istic smallest scales that have boon predicted for differ- 
ent CDM candidates [IO^^Mq for neutralinos (Hofmann, 
Schwarz & Stockcr 2001; Schwarz, Hofmann & Stocker 2002; 
Berezinsky, Dokuchaev & Eroshenko 2003) compared with 
10"^'^ Mo for axions (Kolb & Tkachev 1996)] could have 
observable astrophysical consequences [such as fcmtolensing 
and picolensing of gamma ray bursts (Gould 1992)]. This 
would open up the exciting possibility of discriminating be- 
tween CDM candidates astronomically. 



2 COLLISIONAL DAMPING AND FREE 
STREAMING 

Primordial CDM inhomogeneities on the smallest scales are 
smeared out by coUisional damping and free streaming, if the 
CDM has once been in thermal contact with the hot com- 
ponent of the Universe (as is the case for WIMPs). These 
damping processes give rise to a sharp cut-off in the primor- 
dial CDM power spectrum. In order to make quantitative 
predictions, we assume the CDM particle to be a bino with 
mass m. The bino abundance is fixed at chemical decou- 
pling, which happens at typically Ted ~ m/25, i.e. above 
1 GeV. Below that temperature, elastic scattering of bi- 
nos and fermions (leptons and quarks from the radiation 
component) can still occur and the bino decouples kineti- 
cally at significantly lower temperatures (Schmid, Schwaxz 
& Widerin 1999): Tkd = 10 MeV to 40 MeV, depending on 
the parameters of the supersymmetric model. 

Between chemical and kinetic decoupling the CDM par- 
ticles interact with the cosmic heat bath, which consists of 
all the relativistic particles in the Universe, and this leads 
to coUisional damping of the CDM density perturbations. 
Since CDM is non-relativistic at this epoch, the dominant 
processes are bulk and shear viscosity. Hofmann, Schwarz 
& Stocker (2001) calculated the effect of these processes on 
the primordial density perturbations and found exponential 
damping with a characteristic comoving wavenumber 
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This corresponds to a length scale of ~ 10~^/H at kinetic 
decoupling. The total CDM mass contained in a sphere with 
radius n/kd is Ma ~ 1O~^°M0. 

For a bino, the kinetic decoupling temperature depends 
on m and on the various sfermion masses mj:. Hero we make 
the simplifying assumption that all sfermions have the same 
mass, which is not correct in more realistic models, but is a 
reeisonable assumption for a first estimate. We consider two 
fiducial sets of parameters: model A (B) has m = 100(150) 
GeV and rrif = 230(190) GeV. For these parameters we find, 
neglecting co-annihilations with other supersymmetric par- 
ticles, that chemical decoupling happens at Ted = 4.0(5.8) 
GeV and ilbino ~ 0.31(0.17), thus models A and B have 
CDM densities at the high and low end respectively of the 
range of values found by WMAP. For these models kinetic 
decoupling occurs at Tkd = 33(21) MeV and the damping 
mass scale is Md = 9 x 10"" (6 x 1O"")M0. 

After kinetic decoupling, the CDM particles enter the 
free streaming regime. We calculate the effect of free stream- 
ing by solving the collisionless Boltzmann equation in 
Fourier space on subhorizon scales, to first order in the 
perturbed thermodynamics quantities and neglecting met- 
ric perturbations, starting immediately after kinetic decou- 
pling. We include the spectrum of CDM density perturba- 
tions at this time, and also take into account the matter and 
radiation components of the universe. The net result is again 
exponential damping with a characteristic scale which is in- 
versely proportional to the (time-dependent) free-streaming 
length scale (fefs = 2\/3/Zfs), multiplied by a polynomial 
term which arises from the ratio of the CDM particles ki- 
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netic energy to the thermal averaged kinetic energy. The 
free-streaming scale becomes approximately constant soon 
after matter-radiation equality; 

/ m flcq/akd 
VTkd/ ln(4acq/akd) ao 



-tic. 



1.7 X 10'' (m/100 GeV)^/^ (Tkd/30 MeV)^/^ 



Mpc 



1 -hln (Tkd/30 MeV) /19.2 



(3) 



The damping scale from free streaming depends on cum = 
il^h^ only via the logarithm, we therefore set it equal to 
WMAP's best fit value, ujm = 0.14 (Spergel et al. 2003). 
The corresponding length scale at matter-radiation equality 
~ 10~^ /H and the total mass contained in a sphere with 
radius vr/fcfs is Mfg ~ 1O~^M0. More precisely, for model A 
(B) we find Mts = 9 x 10"''(6 x 1O"'^)M0. 

We can now put together all the damping factors for 
the CDM mass density perturbation A = 5p/p. Well af- 
ter equality (our approximations are valid for a/a^q J? 10), 
when the comoving free streaming length is frozen, density 
inhomogeneities with large wavenumbers are suppressed by 
a damping factor 



Dik) 



(i) " (£) - ii) 



(4) 



The present approximation is valid for fc/fefa < 1, as terms 
of order (fc/fcfg)* have been neglected in the polynomial. 

The mass density contrast as function of redshift and 
wavenumber can now be written 



A{k,z) = A{k,zi)T2^^{k,z)D{k) 



(5) 



where A(fe, Zi) is the primordial density perturbation and 
7a (fe, z) is the transfer function which encodes the gravita- 
tional evolution of A(fc, z). 



3 GRAVITATIONAL GROWTH ON 
SUBHORIZON SCALES 

Following Weinberg (2002), we model the Universe as a 
two component (non-relativistic matter and radiation) fluid. 
Here we utilize the zero shear (or longitudinal or conformal 
Newtonian) gauge. It is useful to define for each fluid (the 
index 'a' stands for radiation, pr = er/3, or non-relativistic 
matter, pcdm ~ 0) the energy density perturbation 

Aa=^, (6) 

(note that Acdm = A) and the velocity perturbations z;a by 
ife 



k 



(7) 



where Ua is the peculiar velocity. 

The Newtonian gravitational potential (j> is defined via 
the metric, which is in the zero shear gauge given by 



ds^ = -{1 + 24)) dt^ + a{l- 24)) SijAx'Ax^ 



(8) 



in the absence of anisotropic stress, where cosmic time is 
denoted by t. Another useful quantity is the hypersurface 
independent quantity C,(= —T^) = A/3 — 4^ (Bardeen 1980, 
1989), which is conserved on superhorizon scales. We use this 
definition to make contact with the WMAP normalization, 
which is given in terms of TZ (Verde et al. 2003). 



Several comments are in order here. First, we work in 
the zero shear gauge, because all subhorizon quantities can 
be interpreted in terms of Newtonian physics, which is not 
the case in the synchronous gauge, where the Newtonian 
gravitational potential is gauged to zero. Second, neutrinos 
are included in the radiation component in order to allow an 
analytic treatment, i.e. their anisotropic stress is neglected. 
This leads to errors of around 10 % (Hu et al. 1995). Further- 
more we ignore the effect of a non-zero cosmological constant 
or curvature (which only effect the evolution of the pertur- 
bations at very late times) and baryon inhomogeneities. At 
early times the baryons are tightly coupled to the radiation 
fluid, and photon diffusion damping rapidly erases small- 
scale perturbations in the baryon fluid at 2: ~ 10^ to 10^. 
On small scales the tight coupling breaks down prior to re- 
combination, and the baryon perturbations grow, however 
Ab >C Acdm still (Yamamoto, Sugiyama and Sato 1997, 
1998). Post decoupling on scales fc > fcb ~ 10"^Mpc~^ the 
residual electrons allow transfer of energy between the pho- 
ton and baryon fluids so that thermal pressure prevents the 
baryon perturbations from growing, until ~ 150 (Ya- 
mamoto, Sugiyama and Sato, 1997; Padmanabhan, 2002). 
As we are interested in CDM perturbations on small scales 
at early times, we can neglect the perturbations in the 
baryon fluid. 

The gravitational evolution of CDM inhomogeneities is 
described by the corresponding transfer functions. For the 
calculation of the transfer functions the equations of motion 
of the related perturbation variables have to be solved. Un- 
fortunately, this is not possible exactly for all times, however 
there are two overlapping regimes for which exact solutions 
exist. The flrst of these regimes is radiation domination, 
Er 3> Ecdm, for which an exact solution can be found that is 
valid for all scales (Schmid, Schwarz & Widerin 1999). In the 
superhorizon limit {k/a ^ H), 4> — > 4>o, Acdm.r — 30o/2, 
''cdm.r —{4'o/2)(k/aH) and C, —?>4)o/2, while in the 
subhorizon limit {k/a 2> H) with x — k/{\^aH): 



Acdm{x) = 

AAx) = 



In a; -I- 7E - - 



Vcdmix) = 



3a/3(^o 



94>o . „, cos a; 

-— cosx, 0(x) = — 300 — ^ 
2 x^ 



(9) 



with 7e denoting Euler's constant. We see that on sub- 
horizon scales during radiation domination the matter per- 
turbation Acdm grows logarithmically while Ar oscillates 
with constant amplitude. Thus for small scales which enter 
the horizon sufficiently long before matter-radiation equal- 
ity the condition CcdmAcdm ^ frAr is satisfied during the 
radiation dominated era. As long as the CDM density per- 
turbations EcdmAcdm dominate the density perturbations of 
the radiation and baryons, the evolution of Acdm is governed 
by (Hu and Sugiyama 1996): 



j/(i+y) 



d^Acd 
di/2 



1 + 



3\ dAcdm 3 

2 



(l-/b)Ac 



0,(10) 



where y = a/uaq and /b — Slb/f^m is the baryon fraction, 
with best fit value from WMAP ft = 0.16 (Spergel et al. 
2003). The exact solution to this equation is a combination 
of Legendre functions of first and second kind, P^{^1 + y) 
and Qi,{y/TTy) with index u{fi) = (^25 - 24/^ - l)/2. We 
smoothly join the small y expansion of these functions with 
the radiation domination subhorizon solution [eq. @] , to ob- 
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tain the normalisation of the solution to eq. llOt . Expanding 
the result for y 2> 1, we find 



+ b 



(11) 



Acdm(j/) = -900 cy"^^ In 

where fc^q = l/(a//)eq, c(//) = r[l + 2u]/{2T''[l + u]) 
and 6(/b) = l/21n(2V3) - 7e - 1/2 - 2/zy - 2r'[v]/r[h'], 
e.g. i/(/b) = 1.80(2), c[i.(/b)] = 1.37(3/2) and 6(/b) = 
-1.57 (-1.74) for /b = 0.16(0). Note that before Zb, CDM 
density perturbations grow as Acdm oc a^^^. Later the 
baryons follow the CDM and the matter fluctations grow 
as a. For the peculiar velocity and the Newtonian gravita- 
tional potential we obtain 



'^cdm(y) 



y d 



2dy 

k ^ 2 



3 / fccq \ / 

4 V y ^ 



l-/b)- 



y 



(12) 



In the following, we omit the subscript cdm. 

For redshifts Zeq > z > Zh (between matter-radiation 
equality and the epoch at which small-scale baryon pertur- 
bations start growing) we find the transfer function for the 
CDM density perturbations for modes which satisfy k > kh 



TA{k,z) = (6c)2 



In- 



V 1+2 / ' 



(13) 



and the transfer function for the Newtonian gravitational 
potential on these scales is given by 



In^ -Kb 



(¥)'(i^)'^"(14) 



The transfer function for the velocity depends on the initial 
time and is therefore not a very useful quantity. 



4 POWER SPECTRA 

In this section we present the dimensionless power spectra 
(defined as Vx{k,z) = [k^ /2n'^){\X{k, z)\^)) normalized to 
the WMAP measurements (Spergel et al. 2003; Verde et al. 
2003). For simplicity we assume that gravitational waves 
have a negligible contribution to the CMB anisotropies and 
that the density perturbations have a Harrison-Zel'dovich 
primordial power spectrum {n — 1). We find for k > kh and 

2eq > Z > 2b 

VAik,z) 



10-M 

Pv{k,z) 
10-M 

X 

10-M 



= 1.06 



D\k) 



In- 



+ b 



(¥)'^'(fc)(^ 



l + z 



In- 



0.60 c 
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+ b 
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(15) 
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where A = 0.9 ±0.1 according to Spergel et al. (2003). Note 
that from the WMAP data n = 0.99 ±0.04, which is consis- 
tent with our assumption of n = 1. The scale of equality is 
feeq = (0.01/Mpc)(a;m/0.14) and 1 + Zeq = 3371(wm/0.14). 

Fig.0shows the power spectrum for the density contrast 
at a redshift of 500, close to the end of the linear regime of 
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Figure 1. The dimensionless power spectrum of the CDM density 
contrast a.t z = 500 for models A and B from the text (full lines). 
Without the effects of coUisional damping and free streaming, the 
power spectra would be given by the dotted lines. 



structure formation. It can be observed that the induced 
cut-off is indeed very sharp and that the power spectrum 
has a maximum close to the cut-off. 



5 DISCUSSION 

The mechanisms of coUisional damping and free streaming 
of WIMPs lead to a cut-off in the CDM power spectrum, 
which sets the typical scale for the first haloes in the hier- 
archical picture of structure formation. A rough estimate of 
the redshift at which typical fluctuations on comoving scale 
R go nonlinear can be made via 



(t(_R, Znl) = 1, 

where (j{R, z) is the mass variance defined by 



a{R,z) = 



W^{kR)VA{k,z)^. 

k 



(18) 



(19) 



where W{kR) is the fourier transform of the window func- 
tion divided by its volume. In accordance with the usual 
procedure, we take the window function to be a top hat. For 
this calculation we need the power spectrum on all scales, 
however our calculation of the CDM transfer function is only 
valid for fc > fcb ~ lO^Mpc"^. We therefore instead use 
the matter transfer function found neglecting the baryon 
density (i.e. /b — 0) which is valid for k > fceq and take 
Jim ~ 0.36(0.22) for model A (B). This introduces errors at 
the 10% level, however the criteria for nonlinearity [eq. I|18|l ] 
is only an order of magnitude estimate. Finally we include 
the power spectrum on scales k ~ fcoq by normalizing <j{R, z) 
to as = (j(8//iMpc,0) = 0.9±0.1 (Spergel et al. 2003), tak- 
ing into account the suppression of the growth of A at late 
times due to the cosmological constant. 

In fig. |2| we plot Znl ; 3'S defined by eq. II18|I . as a func- 
tion of the scale R. The plateau at 7? < 1 pc is due to the 
sharp cut-off in the power spectrum. We can now give a more 
precise picture of the onset of the hierarchical structure for- 
mation process; non-linear structure formation starts at a 
redshift z^{^^, which takes values in the range 30 to 80, and 
2^'"' ~ 60 for the best fit WMAP matter density. To be 
more specific, this is the epoch when the typical overdense 



CDM power spectrum on suh- galactic scales 5 



100 




log,„(R/Mpc) 

Figure 2. The redshift at which typical fluctuations of comov- 
ing scale R become non-linear, for the two models discussed in 
the text. The full lines take into account the effects of coUisional 
damping and free streaming, whereas the dashed lines show the 
behaviour without a cut-off in the power spectrum. 



regions go non-linear. Note, that rare fluctuations with large 
amplitude will go non-linear at much higher redshifts. If the 
density fluctuations have a gaussian probability distribution 
then a Na fluctuation will go non-linear at roughly A^Zni. 
We leave the discussion of these rare large fluctuations and 
their cosmological consequences for a forthcoming paper and 
conclude with some comments on the typical fluctuations. 

Let us estimate the size and mass of the flrst gener- 
ation of subhaloes that form at z^{^^ using the spherical 
collapse model (e.g. Padmanabhan 2002). The mean CDM 
mass within a sphere of comoving radius R is M{R) = 
1.6x 10"'^Mo(a;m/0.14)(_R/pc)^. CDM overdensities that go 
non-linear have mass twice this value i.e. roughly equal to 
the mass of Mars. These WIMP haloes are however much 
less compact than Mars. The physical size of the first haloes 
at turn-around (when their evolution decouples from the 
cosmic expansion) is r = 1.05i?/(l + z";""") ~ 0.02 pc. 
The first haloes then undergo violent relaxation, decreas- 
ing in radius by a factor of two so that their present day 
radius is of order tens of milli-pc (comparable to the size 
of the solar system). If some of these first haloes could sur- 
vive to the present day their overdensity would be of or- 
der Ahaio = 7(1 -I- 2:^™)^ ~ 10*^, several orders of magni- 
tude larger than that of galaxies. Rare fluctuations, a non- 
gaussian fluctuation distribution, or a blue (n > 1) primor- 
dial power spectrum could lead to even larger overdensities. 

We regard this letter as a first step towards an ab-initio 
calculation of the small scale structure of CDM. The result- 
ing power spectra are presented in a form that can be used as 
input for future very high resolution simulations. Only once 
the fate of the first CDM haloes has been understood in 
detail, will it be possible to make robust predictions for the 
expected signals in direct and indirect dark matter searches. 
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